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Processing Seismic Data 

The present invention relates to a method of processing multi-component seismic data. 
It particularly relates to a method of processing seismic data to determine a calibration 
filter that calibrates one component of the seismic data relative to another component of 
the seismic data. The invention further relates to an apparatus for processing seismic 
data. 

Figure 1 is a schematic view of a seismic surveying arrangement. In this figure the . 
surveying arrangement is a marine surveying arrangement in which seismic energy is 
emitted by a seismic source 1 that is suspended within a water column 2 from a towing 
vessel 3. In this example the water column is the sea, but the methods described 
hereinbelow can be applied to data acquired in seawater or in freshwater. When the 
seismic source 1 is actuated seismic energy is emitted downwards and is detected by an 
array of seismic receivers 4 disposed on the seafloor 5. (As used herein the term 
''seabed" denotes the earth's interior, and the term "seafloor" denotes the surface of the 
seabed.) 

Many seismic surveys now use multi-component receivers that record two or more 
components of the seismic energy incident on the receiver. For example a 3-component 
(3-C) seismic receiver contains three orthogonal geophones and so can record the x-, y- 
and z-components of the particle motion at the receiver (the particle motion may be the 
particle displacement, particle velocity or particle acceleration or even, in principle, a 
higher derivative of the particle displacement). In a marine seismic survey a 4- 
component (4-C) seismic receiver can alternatively be used. A 4-C receiver contains a 
(dynamic) pressure sensor such as a hydrophone in addition to three orthogonal 
geophones and so can record pressure fluctuations as acoustic waves propagate in the 
water column (a scalar quantity) in addition to the x-, y- and z-components of the 
particle motion of the seabed. 
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Many different paths exist by which seismic energy may travel from the source 1 to a 
receiver 4 in the seismic surveying arrangement of Figure 1. A number of paths are 
indicated schematically in Figure 1. 

The path 6 shown in Figure 1 is known as the "direct path". Seismic energy that travels 
along the direct path 6 travels from the source 1 to a receiver 4 essentially in a straight 
line without undergoing reflection at any interface. 

Path 7 in Figure 1 is an example of a "water layer multiple path". Seismic energy that - 
follows a water layer multiple path propagates wholly within the water column 2, but 
undergoes one or more reflections at the surface of the water colunui and/or the seafloor 
5 so that the seismic energy passes through the water column more than once. The 
water layer multiple path 7 shown in Figure 1 involves one reflection at the seafloor 5 
and one reflection at the surface of the water column, but many other water layer 
multiple paths exist. 

The path 8 in Figure 1 is an example of a ''critical refraction path". Seismic energy that 
follows the path 8 propagates downwards to the seafloor 5, and penetrates into the 
earth's interior 10 (ie into the seabed). The seismic energy continues propagating 
downwardly, until it reaches a boundary 11 between two layers of the earth that have 
different acoustic impedance. The seismic energy undergoes critical refraction, 
propagates along the boundary 11. before eventually being refracted upv\/^ards towards 
the receiver 4. Critical refraction may also occur at the water-seabed interface, and 
downwardly propagating seismic energy that is refracted in this way will propagate 
along the water-seabed interface and will then propagate upwardly into the water 
coluihn. 

The path 9 shown in Figure 1 is known as a "primar>' reflection path". Seismic energy 
that follows the primary reflection path 9 propagates downwards through the water 
column, is refracted at the seafloor 5, and propagates downwardly through the earth's 
interior. The seismic energy is refracted at the boundary 11, but is not critically 
refracted and so continues to propagate downwardly into the earth. It eventually 
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undergoes renection at a geological stmcture 12 that acts as a partial reflector of seismic 
energy, and the reflected seismic energy is. after further refraction as it passes upwardly 
through the boundarj' 11. incident on the receiver 4. The general intent of a seismic 
survey is to make use of the seismic energy that follows the primary reflection path in 
order to obtain information about the interior structure of the earth. 

Seismic energy acquired at a receiver may contain upwardly and/or downwardly 
propagating seismic energy depending on the location of the receiver and on the event. 
For example seismic energy that travels along the critical refraction path 8 shown in 
Figure 1 will, when it is incident (travelling upwardly) on the water-seabed interface, be 
partly transmitted into the water column and partially reflected back into the seabed 10. 
Thus, a critical refraction event will consist purely of upwardly propagating seismic 
energy above the seafloor 5, but will contain both upwardly and downwardly 
propagating seismic energy below the seafloor 5. As another example, seismic energy 
that travels along the direct path 6 shown in Figure 1 will, when incident on the water- 
seabed interface 5. be partially transmitted into the seabed and partially reflected back 
into the water column. Hence, the direct event will contain both upwardly and 
downwardly propagating seismic energy above the seafloor, but will contain only 
downwardly propagating seismic energy below the seafloor. It is therefore often of 
interest to decompose the seismic data acquired at the receiver 4 into an up-gomg 
constituent and a down-going constituent, above or below the seafloor 5. For example'! 
in a 4-C seismic survey it may be of interest to decompose the pressure and the vertical 
particle velocity recorded at the receiver into their up-going and down-gomg 
constituents above the seafloor. 

Various filters that enable decomposition of seismic data into up-going and down-going 
constituents have been proposed. For example, K.M. Schalkw.jk et al have suggested, 
in "Application of Two-Step Decomposition to Multi-Component Ocean-Bottom Data- 
Theory and Case Study", J. Seism. Expl. Vol. 8 pp26l-278 (1999), that the down-going 
and up-going constituents of the pressure just above the seafloor may be expressed as: 



M&C Folio No P522S7 WO 



where P is the pressure acquired at the receiver. F is the up-going constituent of the 
pressure above the seafloo, is the down-going constituent of the pressure above the 
seafloor./is the frequency, k is the horizontal wavenumber, v, is the vertical panicle 
velocity component acquired at the receiver, p is the density of the water, and g ,s the 
vertical slowness in the water layer. 

As can be seen, the expressions in equation (1) require two of the components of 
setsmic data recorded at the receiver to be combined. These filters are an example 
where it is necessary to combme two components of the acquired sdsmic data. It may 
also be necessary to combme two or more components of the acquired seismic data in 
order to decompose the acquired seismic data into p-wave and s-wave (pressure-wave 
and shear-wave) components, or to remove water level multiple events from the seismic 
data. 

One problem in combining different components of the seismic data acquired at a 
receiver is that the different components of the seismic data may not be correctly 
cahbrated against one another. This is particularly the case wher. the two components 
that are being combined are. as in equation (1). the pressure and the vertical particle 
velocity. Ihere are usually differences in coupling or impulse response between the 
hydrophone used to acquire the pressure and the geophone used to acquire the venical 
parttcle velocity. It is necessary to calibrate the data for these differences before the 
pressure and vertical particle velocity can be combined, and this process is known as 
"P/v. calibration". This calibration process involves developing a calibration filter that 
compensates for the differences in coupling and ii^pulse response between the 
hydrophone and the vertical geophone and then applying the filter to one data set to 
compensate for the differences in coupling. 
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Schalkvvijk et a!, and others, have suggested that the PA, calibration problem can be 
addressed by assuming that one component of the seismic data has been correctly 
recorded, and calibrating the other component of the seismic data against the componen. 
that ,s assumed to be correctly recorded. In general, it .s assumed that the hydrophone 
.s well coupled to the seismic wavefield. so that the pressure recording is taken to be 
con-ect. The vertical component of the particle velocity is then calibrated against the 
pressure to compensate for coupling and impulse response d.fferences between the 
hydrophone and the vertical geophoae. Schalkwijk et al therefore proposed that 
equation (1) above should be modified by applymg a calibration filter to the vertical 
particle velocity. They proposed that the equation given above for the down-goinc 
constituent of the pressure above ihe seafloor should be modified to read as follows: 

In equation (2) a(f) represents a frequency-dependent calibration filter. The remaining 
terms in equation (2) have the same meaning as in equation (1). 

The method proposed by Schalkwijk et al. for deten^ning the calibration filter «(/) ,s to 
find the calibration filter that minimises the energy of the down-going pressure 
constituent above the seafloor for a portion of the seismic data that contains only 
primary reflections. Seismic energy travelling along a primary reflection path is 
propagating upwardly just above the seafloor at the receiver position, so that the down- 
gomg consatuent of the pressure just above the seafloor should be zero for a pnmary 
reflections. Schalkwijk et al. used a least squares method to find the calibration filter 
that mrnimises the energy of the down-going constituent of die pressure in a window 
containing only primary reflection events. Once, the calibration filter a(f) has been 
determined in this way. it is applied to the entire data set to calibrate the vertical particle 
velocity. 

The existence of various paths of seisn.c energy from the source to the receiver means 
that the data acquired at the receiver in a real seismic survey wi„ contain events 
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corresponding lo more than one possible path. These events will occur at different 
times after the actuation of the seismic source 1. as different paths of se.smic cncrov 
have different associated .-avel times. In order to apply the method of Schalkwijk et 7l 
to determine the calibration filter, data in a time window that contams only a primary 
reflection event must be selected. 

The method proposed by Schailavtjk et al. has the disadvantage that the time window 
contammg only primary renection events has to be picked manually. The primary 
reflection events arc not the first events acquired at the receiver following actuation of 
the source, and so cannot be picked automatically. A futther disadvantage is that in 
some cases, for example if the seismic source has a long signature, it may be hard to 
d,sttngu,sh between the dn ect am val and the primary reflection events, so that it may be 
difficult to isolate the correct events. Moreover, in shallow water the water layer 
multiple events may arrive shortly after the direct wave. In this case, the derivation of 
a(0 IS based on a very limited amount of data, reducing the accuracy of the results The 
direct event and water multiple events contain downwardly propagating seisirtic energy 
so that use of a time window that inadvertently included the direct event or water 
multiple events would not give correct results for the calibration filter 

Co-pending UK Patent appUcation No. 0200560.1 and PCT appHcation PCT/GB 
03/00052 propose applying the approach of Schalkwijk et al to a time window 
containing critically refracted waves, which also contain only up-going energy above 
the seafloor. This technique is moi^ suitable for automation using first-break pickers 
and allows the method to be applied to data acquired in shallow waters. 

Ball. V. L. and Cotrigan. D. suggest, in "Dual sensor summation of noisy ocean-bottom 
data". 66th Ann. Intemat. Mtg: Soc. of E.Kpl. Geophys.. 28-31 (J996). calibrating the 
vertical geophone data against pressure by applying a hydrophone ghost operator to v 
and a geophone ghost operator to the pressure. It can be shown that this approach is 
equivalent to predicting the down-going pressure reflected from the sea smface from the 
computed up-going pressure. N^nin.sing the difference between the predicted and the 
computed down-going pressure reflected from the sea surface at times larger than 
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source duration plus the one-way propagation time through the water layer then allows 
the desired calibration filter to be determined. A significant drawback of this techn.que 
is the need for accurate information about the depth of the water column in the survey 



area 



Another approach to calibrating data for differences before the pressure and vertical 
particle velocity ,s to require sdsn.c energy to be preserved during propagation through 
the water layer. In thus case, frequency- and wavenu«.ber-dependent calibraUon 
operators are designed by means of spectral balancing of the up- and down-going wave 
constituents just above the seafloor. 

The present mvention provides a method of processing multi-component seismic data 
acquired at a receiver station from seismic signals propagating in a medium, the method 
compnsing the steps of: selecting a first portion of the seismic data; and dete^ning a 
calibration filter from the first portion of the seismic data, the calibration filter being to 
calibrate a first component of the seismic data relati ve to a second component of the 
seismic data; wherein the step of determining the calibration filter comprises processing 
the data in the common shot domain. 

The tern, "receiver station" may denote, for example, a single receiver (for example in 
single sensor seismic acquisition), a hard-wired group of receivers, etc. 

Coupling and instalment response variations for a particular seismic receiver or receiver 
station are receiver-consistent effects - that is, once a multi-component receiver has 
been deployed the required calibration operator for the receiver is constant and does not 
depend on the source position. Receiver calibration has therefore hitherto been 
perfomred in the common-receiver domain (CRD). When seismic data is processed in 
the common-receiver domain, all data acquired at one receiver is sorted into a "gather" 
for that receiver - and since all data in the receiver gather was acquired at a single 
receiver, a single calibration operator can be used to correct all the data in the gather for 
coupling and instalment response variations. 
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The inventors have realised, however, that optimising the vector fidelity of mult, 
component receiver, in the common shot domain (CSD) has cons.derable advantages 
including: ' 

unlike CRD processing, processing in the CSD does not impliedly assume the 
sub-siirlacc to be laterally invariant; 

rough sea perturbations can properly be accounted for; 

- it is possible to benefit from the dense receiver spacmg if the source-side 

spacing IS coarser; and 

- waterbome no.se that ,s not shot-generated can be removed more efficiently. 

It is convent,onaI to perform the cahbration technique in the common-receiver domain 
smce th. allows the cahbraUon filter for each receiver to be determined separately but 
th. has the disadvantage that rough sea perturbations cannot be properly accounted for 
and that background noise cannot be removed efficiently. Furthermore, applying the 
processmg in the common receiver domain implicitly assumes that the earth's sub- 
surface is laterally invariant. 

The calibration filter may then applied to seismic data recorded at the speafic receiver 
or receiver station which acquired the seismic data used to determine the calibration 
filter. Applying the calibration filter compensates for the effects of the different 
couphng for the first and second components of the seismic data. 

The calibration filter is receiver-specific, so a separate calibration filter rs preferably 
determined for each receiver (or receiver station) in a seismic survey. 

A second aspect of the invention provides a method of seismic surveying comprising 
the steps of: actuafing a source of seismic energy; acquiring seismic data at a receiver 
station spatially separated from the source; and processing the data by a method as 
defined m the first aspect of the invention. 
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A ,h,rd a.p« of invention provides an apparalns for p,^ess,ng n,u,„-c„„,p„„e„, 
scsmic da^ acquired a. a receiver oration fron, .eismic signals propagacn. ,„ » 
™ed,„m, ,he apparain, comprising: means for dctermrning a calibra„„„ f,|,cr from , 
fi... portion of ,he seismic da.a. d. cai.tadon fil,cr bemg ,o calibraie a firs, component 
of the seismic data relative to a second component of the seismic data; wherein the 
apparatus is adapt«l to determine the cahbratton ftlter by processing dte seismic data i„ 
the common shot domain. 



The apparatus may comprise a programmable data 



processor. 



A fourth aspect of the ,„vent.on provides a storage med.um containing a program for 
contromng a programmable data processor to carry out a method as defined in the first 
or second aspect. 

A fifth aspect of the invention provides a program for controlling a computer to carry 
out a method as defined in the first or second aspect. 

Preferred features of the invention are set out in the dependent claims. 

Preferred embodiment, of the present invention will now be descnbed by way of 
Illustrative example with refei^nce to the accompanying figures in which: 

Figure 1 is a schematic illustration of a seismic survey; 

Figure 2 is a block flow diagram of a method accordmg to one embodiment of the 
present invention; 

Figure 3 is a bl«:k flow dtagram of a method aecord.ng to a second etnbodiment of the 

present invention; and 

Figure 4 is a schematic block dtagram of an apparatus according to the present 

invention. 
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A ou.l,„.<, »b„.e, exis.l„g PA, c.,ibra,io„ ,ech„l,„es a. b.ed „„ fi„di„, ca,i..a„„„ 
fters .ha, ™„tai.e ,he „p-g„i„g „, ,„„3,„^„, ^^^^^^^ ^ 

column (a. defined by e<,„«,on (2)). The term i„ 

considered as a spadal fiUer, and ,he pHor methods have the disadvantage ,ha, ,„ev 
no,t„ne computing the te,™ 7/,^ „^ ^^^^^^^^ 

ot pcsttion, spa.,.1 mter operations on the ,e„„ „„p„^^,^ 
computation of the calibration flltet, a(f) for each ,„d,v,dual receiver posiUon. 

Tlte present tnvention ovetcomes tbese dtsadvantages by perfott^ng the spatial filter 
operatton in the comtnon shot domain. This allows the opttmisation criteria to be 
appbed to the up- and dovvn-going components of parade velocity, .bus avotding 
spatial filter operations on the tenn [a(f) vz(f,k:)J. 

The invention makes use of the fact that the up- and down-going constituents of the 
vertical component of particle velocity in the water layer can be expressed as 

2 2p ^ (3) 

In cttation (3). oW again denotes the ft^uency-dependent calibtatton filter that 
con^cts for imperfeedons i„ the reeordtng of v. v, denotes the up-gotng constituent of 
*e verttcal component of pattCe velocity, and v.- denotes the down.going constttuen, 
Of the ve*cal component of pat.Ce vel^ity just above the seafloo. The remaintng 
quanttUes in equation (3) have die same meaning as in equation (2). 

Ue detentrination of the up- and down-going cons„t.,e„,s of the vertical component of 
U^ patttc e velocity involves applytng the spa.,., filter to the pressure, rather than 
10 the caltbrated vertical panicle velocity. As a rcs„„, .he fil.er operation does not 
.nteract wiU. the unknown calibration filter am when apphed rn equatton (3) in the 
common shot dotnain. Detenoimng the calibration filter in .be common sho. domain - 
or perfonning a, least some of .he steps of the de.en„i„ati„n of the calibration filter in 
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.he co„,mo„ sho, domain - allows the calibraUon filter ,„ be deK^ined fron, ,„c „„ 
and down-soing co„,po„e„u of ^ panicle veloei.y according ,o cquaUon (3), fcrcfore 
greatly simplifying the processing. 

A funher advamage of aaing equadon (3) is >ha, ,he filter operator U^f.t, „,ed ,„ 
equalton (2) contains a pole, fn contrast, the filter operator ,m u, et,„atio„ (3, does 
not contain a pole, and sinrply contains a zero. As a resuit, techniques ba^d on 
jatton (3) ate numerically more stable, and fac.litate the use of .spadaliv compact 
fi .er appro^mations. ne use of spattally compact filter approximations is 
advantageous when processing data acquired a. a st^vcy locafon whete thcte are 
significant lateral variations in the earth's surface pat^eters. 

Hgure 2 tliustratcs one method of determining a clibmtion filter in the common shot 
do..,„. This e^ple ,s descnbed with reference to P/v. calibradon, but a similar 
procedure may also be carded out for P/v, calibraUon as suggested by Sehalhvijk et .1 
(above), provided drat the elastic p„>perties of the seabed are known. 

The method shown in figure 2 is based on the P/v. calibration ,^h,u,„es proposed by 
Schalkwijk e, al (above) and in UK patent applicadon No. 0200560.,. I„ the method of 
figu„ 2, the frequency-dependent calibration filter is detennined by mimmtstng the 
eiteigy of the down-going constituent of the verfical component of the panicle veloc.ty 
V. over a data window that contains only up-going events. For example, the data 
vvmdow may com^n only the primary reflecfion even, (as proposed by SchaUrwtjk et 
al.) or ,t may contain only critically refiracted events (as proposed in UK patent 
application No. 0200560. 1 ). 

Inidally. multi<ompo„ent seismic data aie acquited at step 20. In diis embodiment the 
multi-component seismic data contain a, least pressure recordings and veriical p.™icle 

velocity recordings. 

■fhe invention may al.einatively be applied ,„ pre-existing seismic data. In this case 
..ep 20 is replaced by the step, step 21, of retneving multi-component seismic data' 
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containing a, te. pressure recordings and vertical parUele vel«:i,v recordings fton, 

Storage. ' ^ »n 

a. s,ep 22. ,he p,ess„re da,a is sorted in.o con,n,on sho, ga.hers. Each ga,„c 

contain all pressure data records f„r a p„„ieular sho, - so where da,a is ac,„,rcd „s,n, 
.he se,sn.e snrve^ng arrangement shown in f,g„o 1, for example, one ga,her w,h 
contatn pressu„ data teeorded when the source , is actuated at one locat.on ,„d.cated 
schentattcally by ,3 in Kgure ™,her gather w„l conta.n ail pressure data receded 
When the source is at another location .ndicated schematically by 13a and so on Each 
gather wU, contain pressute data ac,„ired at all „ceive. 4 in the setstnie su,.ey,n. 
atrangement. Tbe pressure data ac,u,red by a receiver w,„ typtcally be in the of I 
..cotd of the pressure at arece.ver measured as a f„net,on of ttme since the actuat.on of 
the seismtc source 1. That ts, each rece, ver records P(rece^»r posUon. „ which, if the 
sources and receivers are arranged along a straight line, may be sitttplified to PM 

Next, a, step 23, the pressure data ate transformed, in this embodiment into the 
fre,ue„cy.wave number m domain to give Pa». lite transfonn to the ftequency. 
wave number domatn may be catted ou, using any sui.ble techntque. n,e invention is 
no., however limtted to the/-* domain and the method of the invendon may be used ,„ 

other domains such a., for example, the domain, the x-/dom^n, the dotnain. the 

a:-/ domain, the tau-p domain, etc. 

Next, a, step 24. the ptessute data are spatially filteted. That ,s, the filter ,m is 
app hed to the pressure data to yteld the filter^ pressute data .mPif.., The «,er 
m) .s obtatn^l from the seismic data acquited a, step 20 or retneved a. step 2, and 
thts can be done u,,ing any tatown technique. The de,er„.„at,o„ of the filtered ptessute 
data may also be calculated using any known technique. This spanal filtenng step ,s 
canned out on each ptessure record.ng in a shot gather, and Ms process ,s then repeal 
for all shot gathers. 

A. step 25, the filtered pressure recordtngs are transformed back to the x-domai„ for 
example the f-x, r-x or t-x doma.n. The filtered pressure reeotdings a,e then so,r«l'i„to 
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receiver gathers. Each receiver gather contain.: rfiit^r„^\ 

^1 gdincr contains (tillered) pressured recordings acquired 

at onJy one receiver or receiver station. 

A. step 26, a fcquen.y dependent ealibralion m,er atf) ,s de.crmmed f„, one receiver 
gather from a pomon of .he .eismie d.,a n, ,ha, gather. ,„ ,,,,3 e,nbod,™„, the 
cai,b,..ion filter is determined nsing e,„ati„„ (3), by find,„g the calibration Alter that 
mtntmtses the „p-g„,„g or dow„-g„,„g constituent of the vertical panicle velocity for a 
..n,e Window in which that cons.tuent ts expected to be ..ro. (Although equation (3) is 
wncten in U,e /./- domain, a cotrespondtng equation ex.sts in the x-dom.in ) For 
example, the time window selected may be a time window dta, is expected to contatn 
only up-go,„g seismtc energy just above the sea floor, such as a time window tha, 
eontains only the pnmary reflection events or a time window tha, contains only 
cntic^ly refracted events. Where a dme window that is expected to contain only up- 
sotng seistntc energy is .^elected, the ealibmUon filter Is detennined by finding the filter 
that mtnimtses me down-going eonstifen, of dte venical panicle velocity in the 
selected time window. 

The result of step 26 is a calibradon filter for the pa«icular receiver (or receiver statton) 
used to acquire Ute pressure data in the n=ceiver gad,er on which step 26 was perfonned 
■nns calibraUon filter may then be used to calibrate verfieal velocity component data 
acquired a. that receiver (or t^iver sution), a. step 27. The cahbration filter 
determined for a receiver (or .ceiver stadon, may be apphed to venica, velcty 
cotnponent data outside the time window used ,o determtne the cahbratton filter and it 
may be used to calibrate veraeal velocity component data which was acquned at that 
receiver (or teeeiver station) but wh.ch was not used ■„ the determination of the 
calibration filter. 

steps 26 and 27 may then be repeated for other receiver gather,, ,„ determine 
cahbratton filters for each rece.ver (or receiver station, in the setsmic surveyin. 

arrangement. 



I 
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The calibr.>.d ve„.al ™,oci.y da,a may b= subjcc^d ,0 further prr»e..i„g srcp. („oi 
shown). 

In .h. =mbodl™e„. of figure 2. *p 22 conpnVs >ra„.,fon„i„g .be pressure d,u:, ,„ ,be 
f-k domain. I, would al«™a,ively be possible for step 22 ,„ compose ,ransf„™,n, ,be 
pressure dau .o ,he r-p do™,„, ,„ wbieb ease step 23 would eontprise fil,e„„:,be 
pressure data in the t-p domain. 

A second en,bodin,e„l of ,i,e present ,„ve„tron is tliustrated in figu,. 3 of the pt^nt 
appbeation. This embodiment is based „„ the principle that the water column is non- 
attenuauve for se.smic waves, and that the free sutface of d,e water column has a 
known, constant reflectton coefficient r„. As a resui,, al, up-going energy i„ the water 
column should, at some later t,me, be recorded as down-going energy but with an 
atttphtude reduced by a factor r. The only e^cepbon to this principle is the diteet wave 
6 - the direct wave 6 produces down-going energy i„ the water column .hat does not 
have any cotresponding up-going energy. A receiver-consistent, ftequency-dependent 

calibration niter can be de.em.ued by ftndingthe calibration Alter thatprovides spectra, 
balancing between .he up-going energy and .he dowu-going energy (after removing or 
muung", ,he direct wave front d,e down-going energy). Thus, the calibratton ftltei is 
determined as ttte ealibraUon filter ttta. minitnises fte following objecive fuucion 



£=Zi»'(/,*)(|v.-(/,t))|-)v,>'>-(/.t,)|y, 



(4) 



In equadon (4), Wm is a weighting funcfcn .ha, spccftes .he fre<,„ency.wave number 
»,ndow over which *e tt^nimisarion process ,s performed. V. ,s the up-going 
constituent of the ver^cal particle velocity, and ,s the down-going constituent of 
the vertical particle velocity after muting of the dtrect wave. 

The terms V.m and V'"' ^ ,„ c,uatio„ (4, ,„ay be expanded us.ng equation (3) to 
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The superscript (r) again indicates that the di,«:t wave has been i^moved. 

M„dng the ditec, wave is a taown. stratghtfonvari p^cess. In essence, a data „ask ,s 
eonstntcted that has a value of zero for the app..,in,ate at^val ,„„es of the d.rec, wave 
and that is equal to one at all other times. The di„=c, wave is then tnuted by 
multiplication of this mask with the data. The amval fme of the direct wave can be 
estimated from a rough esh„,a,e of the survey geometry, water depth and the velocity of 
setsmic energy. The transiUon between zero and one ,n the data ntask is preferably 
smooth rather than abr.pt, and can be anywhere between the arrival t,me of the direct 
wave itself and the arrival ,in,e of the first water layer multtple since any prin^ 
reflections that nray occur between these times are purely up-going waves. The direct 
wave is most easily muted in the (x-t, or (tau.p).domains - in which case dre pressure 
data and vertici velocty data are muted and then transfonr^d .„ the/.* domain. 

Figure 3 is a flow char, showing one example of d,is medrod. WHally muld- 
componen. setsmtc data are acquired at step 30 or are retrieved from storage at Itep 31 
The mult-component setsmic data contain at least pressure tecording. ami vertical' 

particle velocity recordings. 

At step 32. the pressure recordings are sorted into common shot gathei.. At step 33 the 
pressure recordings are transformed to the f-k domain, and at step 34 the filter ,(f.k) ,s 
Obtained and the pressure data .e filtered to determine the quantity steps 
30-34 correspond generally to steps 20-24 of the method of figure 2. 

At step 35. the filtered pressure data are sorted into common receiver gathers. 
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one pa.icu,a. rcccver ga.he. .s .hen selected at step 36. The pressure data , 
the ve^cai ve.ocit, data ^or that gather a. then .uted to re.ove the .1 

je. The ..ed p.ss.e data and ve.ea, v.ocit, data a. then t.ns Wd to , ^ 

on.a.n. at step 37, to g.ve and V<^>,,, ^ 

are also transformed to the f-k domain in this step, to gtve K m 

At step 38. a frequency-dependent caiibrat.o„ filter . dete.n.ned fo. .he selected 
-^ver sather. This calibration .ter . dete^ned as the cal.ration f i ^ ^ 
n^mmtses the objective function E of e,aat.on (4) or ecjuat.on (5). 

The calibration filter obtained in step 38 may then be used to calibrate vertical veloc.y 
data acquired by the receiver (or receiver stnrinn^ ^ 
™ . . . ^. , ^ coirespondmg to the selected gather 

This IS mdicated as step 39 in figure 3. 

Steps 36-39 are then zepeated for other i 



for each receiver. 



receiver gathers, to determine calibration filters 



^rfo^ed ,„ way. and i, „„„„ ^ p„.i,. , ^^^^ 

... ... app„„, an, ,^ ^ 

respective receiver. 

1^ .nve„,.„„ .a. ^„ ^^^^^^^^^ 
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Hsuro 4 i= » sche^a^c block diagram of a programmable apparatus 14 accord,„. ,„ rte 
pr«en, inventa. I„e appara»,a comprises a programmable da.a p^cessor , 
programme memory .6, for instance in ,he form of a read-only memory (ROM,, s.orinl 
a progranrme for controlling the data proc^sor 15 to perform any of the pr^xssing" 
methods described above. The apparatus f„„l,er composes „„„.volat,le read/write 
memory 17 for storing, for example, any data wh.ch mus, be ret.tned in the absence of 
power supply. A "working" or scratch pad memory for the data processor is provided 
by a random access memory (RAM, ,S, An i„p„t mterface 19 is prov.ded. for instance 
for receiving commands and data. An output m.erfacc 20 is provided, for instance for 
d,spl.y,„g infon^ation relattng to the progress and result of the method. SeisnUc data 
for processing may be supplied v,a the ,np„, mterface 19, or may alternatively be 
retrieved from a machine-readabJe data store 21. 

The progratttme for operadng the system and for performing the method described 
herembefore is stored in the programme memory .6. which may be embodied as a semi- 
conductor memoty, for instance of the well-known ROM type. However the 
programme may be stored in any other suitable storage medium, such as magnetic'data 
earner 16a. such as a "floppy disk" or CD-ROM 16b. 
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1. A method of processing multi-component seismic data acquired at a reccvcr 
station from .eismic signals propagating in a medium, the method compns.n, ,he steps 
of: selecting a first portion of the seismic data; and determmmg a caUbrauon f„,er from 
the first portion of the seismic data, the cahbrat.on filter be.ng to cal.brate a first 
component of the se.smic data relative to a second component of the seism, data- 
wherem the step of detemnning the calibration filter composes processing the data in 
the common shot domain. 



2. A method as cla,med in clann J and comprising applying the calibration filter to 



the second component of the seismic data. 



3. A method as claimed m any preceding claim wherein the seismic data are 
obtamed from seismic signals propagating in a water column. 

4. A method as claimed in claim 1. 2 or 3 wherem the second cor^ponent of the 
seismic data is pressure. 



5. A med,od as claimed in claim 4 and comprising determming ,(f,lc)PaJc) in the 
common shot domain, where P denotes pressure. , denotes vertical slowness, /denotes 
frequency and k denotes horizontal wavenumber. 

6. A method as claimed in claim 4 or 5 wherein the first component of the sersm.c 
data is a component of the particle motion. 

7. A method as claimed in claim 6 wherem the first component of the seismic data 
is the vertical component of the particle motion. 

8. A method as claimed in claim 3 or ,n any of claims 4-7 when dependent from 
claim 3 wherein the step of deternnning the calibration filter comprising determining a 
calibration filter that n.nm.ises the energy unmediately above the base of the waL 
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column of the down-going constituent of the first component of the seismic data for the 
selected portion of the seismic data. 



9. A method as claimed in claim 8 and comprising determining a calibrat.on filter 
that minimises the energy immediately above the base of the water column of the down- 
going constituent of the vertical component of the particle velocity accordmc 



to 



<(/.^)=-«(/K(/.*)+^P(/,it) 



where v.* denotes the down-going constituent of .he vertical component of the particle 
velocity. V, denotes the acquired vertical component of the particle velocity. a(f) 
denotes the calibration filter, and p denotes the density of the water column. 

10. A method as claimed in claim 3 or in any of claims 4-7 when dependent from 
claim 3 wherein the step of determining the calibration filter comprising detennining a 
caUbration filter that minimises the difference, after muting the direct wave, betwL 
up-going energy in the water column and down-going energy in the water column. 

11. A method as claimed in claim 10 and comprising detennining a calibration filter 
that minimises the following objective function: 

where Wm is a weighting function, v,' is the up-going constituent of the vertical 
particle velocity, and v «^ ^own-going constituent of the vertical particle velodty 
after muting of the direct wave. 

12. A method of seismic surveying composing the steps of: actuaUng a source of 
seismic energy; acquiring seisnuc data at a receiver station spatially separated from the 
•source; and processing the data by a method as defined in any of claims 1 to 1 1. 
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1.3. A method as claimed in any preceding claim and comprising applying the 
calibration filter to a second portion of the seismic data acquired at the receiver station. 

14. An apparatus for processing mulU-component seismic data acquired at a receiver 
station from seismic signals propagating in a medium, the apparatus comprising: means 
for determining a calibration filter from a first portion of the seismic data, the 
calibration filter being to calibrate a first component of the seismic data relative to a 
second component of the seismic data; wherein the apparatus is adapted to determine 
the calibration filter by processing the seismic data in the common shot domain. 

15. An apparatus as claimed in claim 14 and adapted to filter the second component 
of the seismic data in the common shot domain. 



16. An app^u-atus as claimed in claim 15 and adapted to filter the acquired pressure 
seismic data in the common shot domain. 



17. An apparatus as claimed in claim 16 and adapted to determine q(f.k)P(f.k) in the 
common shot domain, where P denotes pressure, q denotes vertical slowness, /denotes 
frequency and k denotes horizontal wavenumber. 

18. An apparatus as claimed in any of claims 14 to 17 and comprising a 
programmable data processor. 



19. A storage medium containing a program for the data processor of an apparatus 
as defined in claim 18. 



20. A storage medium containing a program for controlling a programmable data 
processor to carry out a method as defined in any of claims 1 to 13. 

21. A program for controlling a computer to carry out a method as defined in any of 
claims 1 to 13. 
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Abstract 
Processing Seismic Data 

A method of processing multi -component seismic data acquired at a recei\'er station 
comprises determining a calibration filter that caHbrate a first component of the seismic 
data relative to a second component of the seismic data in order to compensate for 
differences in coupling between the two components. The determination of the 
calibration filter comprises processing the data in the common shot domain. This 
allows optimisation criteria to be applied to the up- and down-going constituents of 
particle velocity, and this is simpler than prior ait methods of applying optimisation 
criteria to the up-going and down-going constituents of the pressure. 



[Figure 2] 
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